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Abstract, d-dimensional Levy processes can be considered as rough paths 15 with the iterated 
integrals understood in a suitable geometric sense; cf. Williams |23j. The expected value of the 
signature (i.e. the ensemble of all iterated integrals over a fixed time interval), which generically 
describes its law, is shown to admit an explicit expression. In fact, the resulting formula is 
precisely of Levy— Kintchine type , with its (d-dimensional) argument omitted. 



1. Introduction 

Consider a d-dimensional Levy process X with triplet {a,b,K). The law of X at time T is an 
infinitely divisible measures on R"* (and conversely any such measure gives rise to a Levy process ) 
whose Fourier transform is given by the celebrated Levy-Kintchine formula (see e.g. [9l [5l [2T1 [T]) 

Ee*(«,x^> ^ exp |r (^-i {u, au) + i {b, u) + J (exp {i {u , y)) - l-^ {u, y) I{|y|<i}) K (dy)^ | . 

If the Levy measure K has exponential moments, an analytic continuation is possible and the 
Levy-Kintchine formula takes the form 

(1.1) Ee<"'^-> = cxp |t Q {u, au) + (6, u) + J (exp ((y, u)) - 1- {y, u) I{iyi<i}) K (dy)^ 

We recall that the signature [resp. expected signature] of a path [resp. process] of suitable 
regularity is the ensemble of all iterated integrals [resp. expectations thereof) and contains signifi- 
cant information of the original object, somewhat similar in spirit to Fourier series, and has gained 
increasing importance in pure and applied mathematics, e.g. [151 [H] and [HIW. For instance, 
in the last reference these ideas were used towards a novel numerical method, popular in financial 
engineering, for second order parabolic equations. Curiously enough, the following result on the ex- 
pected signature of a multidimensional Levy-processes also answers the (at first glance non-sensical) 
question what meaning is left to (11.11) after dropping (!) all u's from its right-hand-side. 

Theorem 1. Assume K has moments of all orders. Then the expected signature 



E 



E 

n>0 



odXt^ (g) • • • (g) odXt. 

0<ti<---<t„<T 



e T((R'*)) 



with values in the infinite- dimensional tensor- algrehra r((R'')) equals precisely the right-hand- side 
of C3P upon dropping all u's in By this we mean that {u,au) is replaced by the (symmetric) 

2-tensor a, {b,u) by the 1-tensor b and so on; all exp{...} functions are then understood via their 
usual power series in T{{'R.'^)). 
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A few comments are in order. 

(i) The simplest way to understand the iterated integrals which appear above is in the "Wong- 
Zakai" sense, i.e. as the mesh-to-zero limit of the iterated integrals in which the (discontinuous) 
sample path X is replaced by a piecewise linear approximation based on a finite partition of [0, T]. 
This type of integration is also called of Stratonovich type (e.g. in [13] which deals with general 
semi- martingales) or in the sense of Marcus, [1] Ch.6], where a more intrinsic definition is given: 
in essence, time is strechted at each jump which is then replaced by a straight line; in the end, the 
size of this OstretchedO time interval is taken to zero. This point of view was also taken in [23] . 

(ii) A geometric p-rough path, p G (2,3), in the sense of Lyons [15] associated to such a Levy 
process X was constructed in [53]; the iterated integrals can then be understood by a basic result 
in rough path theory (e.g. [15] First Theorem] or [TT] Ch.9]. Continuity of this lifting procedure in 
rough path metrics will be useful in the sequel. 

(iii) Theorem[l]is far from any sort of corollary of the usual Levy-Kintchine formula. To appreci- 
ate the content of this type of result consider first K — Q (case essentially contained in [3 HI]), so that 
X can be written in terms of a d-dimensional standard Brownian motion _B, that is, Xt = ht + aBt 
with (Tcr^ = a. Then 

(1.2) EXt = 6T, 

i-T j,2 i-T 

E/ Xt®odXt = —{b'S>b)+'E aBt®od{aB)^ 
Jo 2 Jo 



{b®h) + - {aB.,aB.)j. 



rp2 I 

— {h®b) + -a 



which is precisely consistent with 



exp |r Qa + 6^ I = exp |t ^ ^ a'^^e, ® + ^ b^eu 



n>0 ' y i,j=l k=l J 

projected to tensor-level 1 and 2. Clearly, any direct computation of the expectation of n-fold 
iterated integrals is prohibitive and underlying algebraic relations between iterated integrals, in the 
spirit of Chen [6] [7] must be used in a crucial way. 

(iv) Its own interest aside. Theorem [1] is effectively the starting point for Cuhature for Levy- 
driven SDEs , in the spirit of Lyons- Victoir [T^], but now with jumps. (Needless to say, jumps are 
widely needed in application areas, notably in quantitative finance.) A curious feature then is that 
the resulting weak approximation is effectuated by (typically: piecewise smooth) paths without (!) 
jumps. Explicit computations of such cubature formulae are possible in many examples; this will 
be discussed systematically elsewhere. 

(v) The afore-mentioned cubature method relies on the fact that the (truncated) expected sig- 
nature arises naturally when one expands the expected value of a time-T-functional of a diffusion 
processes ("stochastic Taylor expansion"). This already makes it clear that the expected signature 
contains a significant amount of information about (the law of) the driving signal. But much more 
can be said: thanks to the shuffle-product [20], [ITl p. 50] the expected signature also provides all 
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moments, say of the finite-dimensional projections, of the random signature. Call /i the law of such 
a projection; if the moments do not grow too quickly - and in particular if the support of /i is 
compact - they will determine /i and then the law of the signature. (A point of view repeatedly 
advertized by Lyons, see e.g. [8] and also [16].) In our situation, exponential integrability of K will 
entail that the law of any such projection enjoyes exponential integrability, and thus (see [4] and 
the refences therein) yields a sufficient condition under which the expected signature described the 
law of the signature. Note that it is a different question, if (or when) the signature [resp. the law of 
some random signature] determines the path [resp. the law of the process]; see [T21[T3] for results 
in this regard. But as already pointed out, as far as many applications are concerned, notably for 
numerics for SDEs, the expected signature contains all the relevant information. 

2. Expected signature of Levy processes 

We start with a (new) short proof of the following classical result (due to Fawcett 0, see also 
[19], [2 Sec.l]). Our argument emphasizes the role of IID increments of Brownian motion (and 
thereby naturally points to Levy processes). 

Theorem 2. Assume X := aB is d- dimensional Brownian with covariance matrix a = crcr-^; 
equivalently X is Levy with triplet (a, 0, 0). Write S {B)^ j, for the sum of all iterated Stratonovich 
integrals, viewed as T (^(W^^^ -valued random variable. Then 

Proof. Set (j)^ := ES {X)^ ^. By Chen's relation [6] and independence of Brownian increments, also 
valid for the stochastic integrals built on these increments. 

Since (t)t ® (t)^ = (t^s ® 4^1 have [0^, = 0, 

log(/)t+^ =l0g(/)t+l0g</)3. 
For integers m,n we have \og4>„i = ^^'^z4>m/n and log0„j — m log It follows that 

log 04 = Hog 01; 

first t = G Q, then any real t by continuity. On the other hand, for i > 0, Brownian scaling 
implies 

where 5\ is the dilation operator, which acts by multiplication A" on the n*'' tensor level, (R'') 
It follows that 

log0i6(R'')^' 

and it remains to identify log^j^ = 1^2 (log0i). To this end it suffies to compute the expected 
signature up to level two; indeed we have already seen in the introduction, equation (II. 2p with 
& = 0, that 

E ("l + Xo,i + / X® odx\ = 1 + ^ X! "^"■''^^ 







. J*'"* 6i >^ 07 . 



Taking the logarithm (in the tensor algebra truncated beyond level 2) then immediately gives the 
desired identification. □ 
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Theorem 3 (Expected signature, Levy no jumps). Assume X is d- dimensional Levy with triplet 
(a, 6,0). Then 



Proof. We offer three different arguments. 

First argument: The previous proof of the drift-free case applies until the conclusion that log 
ES {X)q rp is linear in T . But the remainder of the above argument, based on Brownian scaling, 
fails when b ^ Q. We have Xt = ht -\- aBt where aa^ = a where S is a standard Brownian motion 
under some measure P, say. Then 

E'[5(X)o,t] =exp(CT) 

for some C € r((R'')). We will show that Tr^ (C) = for all fc > 3, the result then follows from 
looking explicitly at levels fc = 0, 1, 2; sec (|1.2p .To this end, we use 



t^o t 



TTfc (C) = lim 



E^[g^] 



t-i-O t 



where g^ := tt^ (S'(X)o,t) for fc e {0, 1, 2, . . . }. Now by Cameron-Martin's theorem, (Xt : t e [0, T]) 
is a Brownian motion under a new probability measure Q. From the scaling of Brownian motion 
as before, we have that for all fc > 3 



lim ■ 



9t 



t2 



= 



But for any t e [0,T], 



E^[g^] = E^ 



and so by Cauchy-Schwarz inequality, 

|E'[5f]|<E^ 



E*- 



9t 



Of course, dV/dQ (which is log-normal) has finite second moments with respect Q and we easily 
conclude that tt^ (C) — for all fc > 3. 

Second argument, along [2,. (Sketch) Use the fact that the signature satisfies a (Stratonovich) 
SDE on the full tensor-algebra. Rewrite in Ito-form with Ito-Stratonovich correction, upon taking 
expectations, one ends up with a linear ODE on the tensor-algebra. It is explicitly solved and gives 
the desired result. 

Third argument: Consider the deterministic "drift" process D : t ^ ht which is Levy with 
triplet (0, 6, 0). It is elementary to check that 



ES {D)q j. = exp 



V fc=i / 
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Then Xt = crBt + Dt = crBt + ht which is Levy with triplet (a, 6,0). The signature Z :— S {X) 
satisfies a hnear Stratonovich SDE of the form 

dZ = Z® odX 

= (z) o dx' 

= V^ (Z) o d {aBf + W^ (Z) b'dt 

with vector fields Vi (z) := Wi (z) :— z ® ei (and omitted throughout). With the usual iden- 
tification of Stratonovich SDEs as RDEs driven by (Stratonovich) enhanced Brownian motion, we 
can now resort to splitting results for RDEs [10, to see that 

S{X)q j, = Zo,T = limS'((TS)o (g)S'(D)o ,r/„ (g) ... 5 (crB)^^ ® (£1)^2^ 

We may exchange limit and expectation (justified e.g. by a Wiener-Ito argument: after all the 
claim convergence is for all tensor-levels, each of which takes places in some inhomogenous Wiener- 
Ito chaos. As is well known, on some fixed chaos, convergence in probability is equivalent to 
i'^-convergence, any q < oo.). Then 



ES (X)o.T 



lim exp 




(gin 



exp 



k=l 



exp <^ T 



2 ^ 



fc=i 



where we used the Trotter product formula in the last step. 



□ 



Theorem 4 (Expected signature, pure jump finite activity Levy). Let X be d- dimensional pure 
jump Levy with triplet (0,0, if). Assume K has finite total mass m — K (R''), in other words X 
is compound Poisson. Assume furthermore that K has moments of all orders. Then 



E5 {X\^ = exp |r j (exp {y) ^ 1) K (dy)| 



Proof. Xt — Ji + ... + JNt is compound Poisson where {Nt) is Poission with rate X = K (R''), the 
Ji's are IID distributed according to K/X. The basic remark is that conditional on n jumps over 
[0,1], 



E 



S{X)^^,\N,= 



and then 



E 



Six)., 



Ve-^^Ei 

n>0 



?i>0 



= e"^ exp (AEe-'i) = exp (A (Ee-'i - l)) . 

Here we used that scalar multiplication with e^^ is consistent with tensor multiplication of e""^-*-, 
an element in T((R'^)) which commutes with all others. Remark also that 



Ee-'i = 1+EJi + -EJl 

so that the expected signature is (component-wise) well-defined provided J (equivalently: K) has 
moments of all orders. □ 



1. 



®2 
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Theorem 5. Let X be d- dimensional Levy with triplet {a, b, K), K finite mass and moments of all 
orders (as before). Then 

ES = exp |t |] a^'^ei ^^j+Yl ^'^'^ + / ("^^P (2^) '^^^ (^2/) j | • 

Proof. Call Y the pure jump components, B a Brownian motion with drift so that X = B + Y, 
note that B and Y are independent. As before, the signature satisfies a linear Stratonovich type 
SDE of the form 

dS {X) = S{X)(E) odX 

= S{X)® {odB + dY) , 

and it follows from splitting 

S (^)o,T = lim-^ (^)o,r/„ ® ^ (^)o,t/„ ® - ® ^ (S)ti^,t » S {Y)^^^^^ 
Exchanging limit and expectation then leads to 

d d 



ES{X)q j, = lim ( expj ^ ( ^ X! a*'^ei®ej +^6*^6^1 | Oexpj^y^ (exp(y) - 1) K (dy)^ 



and the result follows again by Trotter's product formula. 



□ 



Remark 1. A detailed analysis based on conditioning of the jumps seems difficult: condition on 
# jumps and jump sizes and times in [0,T]. Then have diffusion over [0,Ti), a jump of size J\, 
diffusion over [Ti,T2) etc etc 



E 



^ {■^)o,T Wt = n, Ti, . . . ,T„; Ji, . . . , J„ 



E 



5(X)o,i.|A^T = n,Ti,...,T„ 



e*^^! (8)E(e'^i) (X) e 



C(T2-Ti) 



iE(e^")®e^(^-^") 



Under the conditioning Nt = 1, the Ti...T„ are uniformly distributed on the simplex A„ (0,T) = 
{0 < ti < . . . t„ < T}. And hence 



I 

J{0<ti< 



E S(X)r,rr,\NT = n = I dti...dt„e'^*i«)E(e-^i)®e^(*^-*i)(8)...(8)E(e"'")®e^(^-*"l 

<...t„<T} 

Write E (e'^) := E {e'^'), independent of i. By the Campell-Baker-Hausdorff formula 

e*^*! (g) E (e-^i) (g) e'^(*="*i^ ... E (e"^") (g e'^^^"*") 
= exp (CT + n log E (e-^) +...) 

where ... contains iterated Lie brackets o/logE(e"^) and C. From here on it seems difficult to 
continue explicit computations. 

Theorem 6. The previous result remains valid when the assumption K (R) < oo Z5 replaced by 
the assumption that K integrates \y\ near 0. It also remain valid, subject to including another 
compensator term ("exp{y) — 1— I{|j^|<i}y when K integrates \yf' nearO. 
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Proof. Consider a general Levy measure, with finite moments of all order, say K, and consider 
the Levy triplet {a,b,K). Truncation of small jumps gives a family of finite Levy measures, say 
(if : e > 0) given by 

(dy) :=if(dy) xl{|^|>,}. 
For fixed e > 0, consider with Levy triplet (a, b'^, K'^) where 



h\v\<i}VK' {dy) . 



By our previous result, 
ES" y = exp i T 



a^'^e, ®e,+Y^ b'^ek + / (exp {y) - 1-I{|y|<i}y) (dy) \ 
k=i J ) 



As e — )> 0, we have 

J (exp(y) - 1-I{|j,|<i}y) if^ {dy) J {exp {y) - l-^y\<i}y) K {dy) 
which settles passage to the limit on the right-hand-side above. On the other hand, 

[s{X%^:e>0) 

is (for every tensor-level ...) Cauchy in probability; this is a consequence of convergence of {X'^ , 'K^) 
converges to a limiting rough path {X, X) in p- variation rough path sense, any p > 2; cf. |23j . Note 
that all X = J X ig) odX (similar X*^) are understood in the "geometric" sense of Marcus. This also 
clarifies the meaning of 5 {X), similar S {X^), which is understood as iterated integral in the sense 
of Marcus. It remains to check uniform integrability (for every tensor level) to have convergence. 
To this end, we show uniform integrability, in e £ (0, 1], of each component of 5* {X"^), i.e. of 



{0<ti<...t„<T} 



odXf 1 • • • o dX[' 



where I — (ii, . . . , i„) is a fixed multi-index. In fact, we shall establish uniform L^-integrability. To 
this end, set J .= I and note 



E 



re,K 



K shuffle of /,./ 

where the sum is taken over finitely many "shuffles"; see [T71 p. 36]. But then 



Eix: 



e,I |2 

o.tI 



T 



\K\<2' 
In 



< 



J2 ^fcexp|TQa + 6 + J {exp {y) ~ l-^yi<i}y) {dy)^ 



T<")(R'i) 



and it is easy to see since {dy) K {dy) x I{|y|<e} and K integrates y^ near the origin that this 
estimate can be made uniform in e > 0. It now follows that 



ES* (X^)o J. ^ F,S {X)g j, as £ i 0. 
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Conclusion is, for the expect signature of a general Levy process, subject to existence of all moments 
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